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m . 

O . A collective treatment of the 7 = scalar mesons below 2 GeV [cr (550), /o(980), / (1370), /o(1500) 

and /o(1710)] in a non-linear chiral Lagrangian framework that is constrained by the mass and the 
£SJ ' partial decay widths of the / = 1/2, 1 scalars [«;(900), Kg (1430), ao(980) and ao(1450)] is presented. 

The sub-structure of these states in terms of two and four quark components, as well as a glueball 
component is explored, and its correlation with the mass of /o(1370) is studied. Consistency with 
the available experimental data suggests that the a (550) is dominantly a non-strange four-quark 
state, whereas the sub-structure of other 1 = states are sensitive to the input mass of /o(1370). 
C ■ This investigation estimates the scalar glueball mass in the range 1.47-1.64 GeV. 
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CN ' PACS numbers: 13.75.Lb, 11.15.Pg, ll.80.Et, 12.39.Fe 

CO 

I. INTRODUCTION 

CN . 

Scalar mesons play important roles in low-energy QCD, and are at the focus of many theoretical and experimental 
investigations. Scalars are important from the theoretical point of view because they are Higgs bosons of QCD and 
induce chiral symmetry breaking, and therefore, are probes of the QCD vacuum. Scalars are also important from a 
■ phenomenological point of view, as they are very important intermediate states in Goldstone boson interactions away 
Q_i from threshold, where chiral perturbation theory is not applicable. There are 9 candidates for the lowest-lying scalar 
mesons (m < 1 GeV): /o(980) [I = 0] and ao(980) [I = 1] which are well established experimentally 0; cr(560) or 
/o(600) [I = 0] with uncertain mass and decay width p|; and k(90CT) [I = 1/2] which is not listed but mentioned in 
PDG lj. The k(900) is observed in some theoretical models Q>IEQ> as well as in some experimental investigations 
It is known that a simple qq picture does not explain the properties of these mesons. Different theoretical models 
that go beyond a simple qq picture have been developed, includin g: M IT bag model K K molecule , unitarized 
^ • quark model d Q , QCD sum-rules , and chiral Lagrangians 0,13 13 111 13 13 EE 13 113 

The next-to-lowest scalars (1 GeV < m < 2 GeV) are: ^(1430) [I = 1/2]; a (1450) [I = 1]; /o(1370), / (1500), 
/o(1710) [I = 0], and are all listed in p|. The /o(1500) is believed to contain a large glue component and therefore 
a good candidate for the lowest scalar glueball state. These states, are generally believed to be closer to qq objects; 
however, some of their properties cannot be explained based on a pure qq structure. 

Chiral Lagrangians, provide a powerful framework for studying the lowest and the next-to- lowest scalar states 
probed in different Goldstone boson interactions (7T7t, ttK , Trry,...) away from threshold [3. [Tol. ITTL [T3. [l3L fT5| . In this 
approach, a description of the scattering amplitudes which are, to a good approximation, both crossing symmetric and 
unitary is possible. To construct scattering amplitudes, all contributing intermediate resonances up to the energy of 
interest are considered, and only tree diagrams (motivated by large N c approximation) are taken into account. In this 
way, crossing symmetry is satisfied, but the constructed amplitudes should be regularized. Regularization procedure 
in turn unitarizes the scattering amplitude. By fitting the resulting scattering amplitude to experimental data, the 
unknown physical properties (mass, decay width, ...) of the light scalar mesons can be extracted. It is shown in [T(| 
that there is a need for a a meson with a mass around 550 MeV in order to describe the experimental data on tttt 
scattering amplitude. Similarly, in the need for a k meson with a mass around 900 MeV for describing the available 
data on ttK scattering amplitude is presented. Motivated by the evidence for a a and a k, and taking into account the 



1 GeV) into a lowest-lying scalar meson nonet is investigated in [ll| . In this approach, the non-linear chiral Lagrangian 



experimentally well-established scalars, the /o(980) and the an (980), a possible classification of these states (all below 

itigated in [ll| • 

is expressed in terms of this nonet, and it is shown pT| that the consistency with several low energy processes (7T7t 
and ttK scatterings, and F[/o(980) — > tttt]) requires a scalar mixing angle which is more consistent with a four-quark 
assignment for the lowest-lying scalar states. This model also well describes the experimental measurements for the 
77' — > rinir decay 0, and estimates the total decay width of a (980) to be around 70 MeV consistent with a recent 
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experimental confirmation |18|. This model has also been employed to describe the radiative decays (f> — ► 77777 an d 
cj> — > 7T7r7 in ref. |l9j| . 

In order to study properties of the oq(1450) and the ^(1430), in the framework of a non-linear chiral, a next- 
to-lowest lying scalar nonet is introduced in |l4l. Similar two-nonet treatment of the scalar mesons in the context 
of linear sigma model are investigated in |l5l l2d||. It is shown in [2] that a mixing between the lowest nonet (a 
four-quark nonet) and the next-to-lowest nonet (a two-quark nonet) provides a description of the mass spectrum and 
the partial decay widths of the J = 1/2 and J = 1 scalars [«(900), a (980), ^(1430) and o (1450)]. The case of I = 
states is not studied in jlj| and is the goal of the present work. Specifically, we would like to see if the 1 = scalar 
states below 2 GeV can be described within a framework that is constrained by the mixing scenario of ref. ^| , and 
in addition includes terms that are relevant to I = states only, such as mixing with a scalar glueball. The mixing of 
scalar glueballs with qq scalar mesons has been studied in the literature |2lll22|. However, the importance of mixing 
between the qq and qqqq nonets in ref. 14] suggests that the mixing of glueballs with the four-quark nonet should 
also be examined. Therefore, in the present work we investigate a description of the 1 = states in terms of two 
and four quark components, and a glueball component. Moreover, we investigate the sensitivity of the results on the 
input data. 

After a brief review of the mixing mechanism of ref. [LH ] in Sec. II, we study the mass spectrum and the sub- 
structure of the isosinglet states in Sec. III. A description of the two-pseudoscalar decay widths of these states is 
investigated in Sec. IV, followed by a summary and conclusion in Sec. V. The basic formulas are listed in the 
appendixes. 



II. MIXING MECHANISM FOR 1=1/2 AND 1=1 SCALAR STATES 

In ref. [14( the properties of the I = 1/2 and 1=1 scalar mesons, k(900), Kq (1430), a (980) and a (1450), in 
a non-linear chiral Lagrangian framework is studied in detail. In this approach, a qqqq nonet TV mixes with a qq 
nonet N' and provides a description of the mass spectrum and decay widths of these scalars. This mixing provides 
an explanation for some unexpected properties of the Kq (1430) and the ao(1450), which are generally believed to be 
good candidates for a qq nonet but some of their properties do not quite follow this scenario. For example, in a 
qq nonet, isotriplet is expected to be lighter than the isodoublet, but for these two states pj: 

m[a (1450)] = 1474 ± 19 MeV > m [Kq (1430)] = 1412 ± 6 MeV (1) 

Also their decay ratios given in PDG [l| do not closely follow a pattern expected from an SU(3) symmetry (given in 
parenthesis): 

-0.92 ±0.12 (1.51), — ±f ± = 0.88 ±0.23 (0.55), — ^ ^=0.35 ±0.16 (0.16) (2) 



T [Kq* —> irK] r [a — > 71-77] r [a -> 7T77 

These properties of the -KTq(1430) and the ao(1450) are naturally explained by the mixing mechanism of ref. The 
general mass terms contributing to the 7=1/2 and the 1=1 states can be written as 

£^=1/2,1 = —aTi(NN) - bTr(NNM) - a'Ti(N'N') - b'Tr(N'N'M) (3) 

where M. = diag(l, 1, x) with x being the ratio of the strange to non-strange quark masses, and 0,6, a 1 and b' are 
unknown parameters fixed by the "bare" masses: 



m 2 [a ] = 2(a + b) m 2 [a' ] = 2(a' + b') 

m 2 [K ] = 2a+(l + x)b m 2 [K ] = 2a' + (1 + x)b' 



2rV.i — o^, _i_ ti _i_ ^2\'irn — o„' _i_ /1 _i_ „\h> (^) 



where the subscript "0" denotes the "bare" states (i.e. before the mixing between N and AT' is taken into account). 
Therefore, as N is a four-quark nonet and N' a two-quark nonet, we expect: 

m 2 [K ] < m 2 [a Q } < m 2 ^} < m 2 [K ] (5) 

Introducing a simple mixing 

<^ A1 — 7Tr(AW') (6) 

it is shown in 0] that for 0.51 < 7 < 0.62 GeV 2 , it is possible to recover the physical masses such that the "bare" 
masses have the expected ordering of ©. Therefore in this mechanism, the "bare" isotriplet states split more than 
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the isodoublets, and consequently, the physical isovector state ao(1450) becomes heavier than the isodoublet state 
ifp(1430) as observed in QJ. The light isovector and isodoublet states are the ao(980) and the k(900). With the 
physical masses m[a (980)] = 0.9835 GeV, m[«(900)] = 0.875 GeV, m[a (1450)] = 1.455 GeV and m[JTo(1430)] = 
1.435 GeV, the best values of 7 and the "bare" masses are found in 

m ao = ma , Q = 1.24 GeV, = 1.06 GeV, m K > = 1.31 GeV, 7 = 0.58 GeV 2 (7) 

The decay ratios © are also investigated in [l4( in which the scalar-pseudoscalar-pseudoscalar interaction part 
relevant to the 1—1/2 and 1 = 1 scalar states is given as 

^Ll 172 ' 1 = Ae ahc t def N%<t>ld^{ + CTr(Nd^)Ti(d^) + A' 'e abc e def N /d a d^ e b d^ c + C'Tr(V'9^)Tr(9^) (8) 

where A 7 C, A' and C are unknown parameters fixed by decay properties of the scalars, and <f> is the conventional 
pseudoscalar meson nonet. It is shown in [T^j that with parameters 

A = 1.19 ± 0.16 GeV -1 , A' = —3.37 ± 0.16 GeV -1 , C = 1.05 ± 0.49 GeV -1 , C = -6.87 ± 0.50 GeV" 1 , (9) 

a reasonable agreement on the decay ratios of the 1^(1430) and the ao(1450) consistent with (|2J), as well as the 
expected decay widths of the ao(980) and the k(900) can be obtained. In next section, we include the Lagrangians 
© and JSJl together with parameters Q and 10 as part of the mass and interaction Lagrangian of the 1 = scalar 
states. 



III. ISOSINGLET STATES 



The general mass terms for nonets ./V and N', and a scalar glueball G can be written as: 

tfnass = ^aJs' 1 ~ cTt(N)Tt(N) - dTr(N)Tr(NM) - c'Tt(N')Tt(N') - d'Tr(N')Tr(N'M) - gG 2 (10) 

The a priori unknown parameters c and d induce "internal" mixing between the two 1 = flavor combinations 
[(Nl + iVf )/\/2 and Vf] of nonet N. Similarly, c' and d! play the same role in nonet N'. Parameters c, d, c' and d! 
do not contribute to the mass spectrum of the 1=1/2 and 1=1 states. The last term represents the glueball mass 

term. The term dmals' X is imported from Eq. J3J together with its parameters from Eq. Q. 

The mixing between N and N', and the mixing of these two nonets with the scalar glueball G can be written as 

7=1/2,1 _ pTr^Xj.^') 



C 1 





mix 



eGTr (N) - fGTi (N') 



(11) 

1/2, 1 mixing, 



where the first term is given in © with 7 from 0. The second term does not contribute to the / 
and in special limit of p — > —7: 

- 7 Tr(VV) - pTr(jV)Tr(JV') = 7 e Qbc e ade <V'" (12) 

This particular mixing is more consistent with the OZI rule than the individual 7 and p terms and is studied in |22| . 
Here we do not restrict the mixing to this particular combination, and instead, examine a range of p values. Terms 
with unknown couplings e and / describe mixing with the scalar glueball G. As a result, the five isosinglets below 2 
GeV, become a mixture of five different flavor combinations, and their masses can be organized as 

¥■ 



c 1 





mass 



pl=0 _ 



-^F M 2 Fo 



■-FMi F 



(13) 



with 



/ 



F n = 
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(14) 



where the superscript NS and S respectively represent the non-strange and strange combinations. F contains the 
physical fields 



/ ^(550) \ 
/o(980) 1 
/o(1370) 
/o(1500) 
V /o(1710) J 



= K^Fn 



(15) 
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where K 1 is the transformation matrix. The mass squared matrix is 



M 



V2\2c + (1 - 
7 + P 
V2p 

e 



2{c + dx) V2[2c+ (1 + x)d] 
hx)d] m 2 ao +A(c + d) 

V2p 
1 + 2p 

V2e 



2m\, 



7 + P V%P 
y/2p 7 + 2p 

ml, + 2(c' ± d'x) V2[2c' + (1 + x)d'] 



V2[2c' + (1 + x)d'\ 
f 



f 4(c' 
V2/ 



d') 



e 
V2e 
f 

2ff 



(16) 



in which the value of the unmixed / = 1/2, 1 masses, and the mixing parameter 7 are substituted in from (JJJ). We 
search for the unknown parameters c, c', d, d' , e, /, <? and p in 1|16|) by fitting its eigenvalues to the mass of the physical 
states. We take m a = 550 ± 50 MeV expected from chiral Lagrangian treatment of the irir scattering in [lOj, as well 
as the following experimental values from PDG : 



m[/ (980)] 
m[/ (1370)] 
m[/ (1500)] 
m[/ (1710)] 



980 ± 10 MeV 
1200 1500 MeV 
1507 ±5 MeV 
1713 ±6 MeV 



(17) 



The largest experimental uncertainty is on the mass of /o(1370), and for our initial study we take its central value of 
1350 MeV. The sensitivity of the results on the mass of /o(1370) turns out to be the main source of error and will be 
discussed later. Figure ^ shows the dependency of x 2 on parameter p. We see that for — 0.4GeV 2 < p < the \ 2 is 
very small, but significantly increases outside this interval. For five values of p the result of fits are given in table [I] 
in which the sensitivity of the fitted parameters on p can be seen. As far as fitting the isosinglet masses is concerned, 
all chosen values of p in table [Ogive more or less the same description. Parameters d and d' induce SU(3) symmetry 
breaking and are at least an order of magnitude smaller than c and c'. The fitted parameters determine the rotation 
matrix (|15fl which in turn probes the quark substructure of the scalars. For p = the rotation matrix is given in i|18|) 
in which the errors reflect the experimental uncertainties in l|17|) . For other values of p in table [H the corresponding 
rotation matrices are given in Appendix A. The overall results show that the sub-structure of <r(550) is not sensitive 
to p and is dominantly a non-strange four quark combination udud, consistent with the investigation of ref. |llj| . 
The sub-structure of the /o(980) is more sensitive to p; for p < — 0.2GeV 2 it has a dominant ss component, whereas 
for p > -0.2 GeV 2 the non-strange (uu + dd)j\[2 component dominates. The later case has some support in QCD 
sum-rules ||. The /o(1370) has substantial two-quark components with some glueball admixture, and the /o(1500) 
and /o(1710) have large glueball components. For each value of p, the corresponding glueball mass is also given in 
table[I] showing a variation in the range 1.47-1.60 GeV. As p — > —7 the result indicates that the glueball component 
of the /o(1500) significantly increases, but in this limit the \ 2 is relatively large, and therefore, the result is not as 
accurate as the cases given in table [1] For p = —0.5 GeV 2 (which is close to the limit p — > —7) the rotation matrix is 
given in Eq. (|A1|) . and is consistent with the result of ref. p^ . 



Fitted Parameters p = -0.5GeV 2 p = -0.4GeV 2 p = -0.2GeV 2 p = p = 0.075GeV 2 



c (GeV 2 ) 


0.141 


0.181 


0.281 


0.157 ±0.009 


0.201 


d (GeV 2 ) 


-0.00494 


-0.0106 


-0.0175 


-0.00779 ± 0.00187 


-0.00840 


c'(GeV 2 ) 


-0.0389 


-0.0475 


-0.0356 


-0.0128 ± 0.0052 


-0.00184 


d' (GeV 2 ) 


-0.00138 


-0.00407 


-0.00452 


-0.00927 ±0.00116 


-0.00825 


e (GeV 2 ) 


-0.169 


-0.165 


-0.242 


-0.323 ±0.018 


-0.248 


/ (GeV 2 ) 


0.0941 


0.234 


0.323 


0.115 ± 0.0568 


0.0874 


9 (GeV 2 ) 


1.211 


1.249 


1.080 


1.272 ±0.021 


1.159 


m G (GeV) 


1.556 


1.581 


1.470 


1.595 ±0.013 


1.523 


x 2 


0.0740 


« 


« 


w 


0.0932 



TABLE I: Best numerical values for the unknown parameters in Lagrangian ifTSy with m[/ (1370)] =1.35 GeV and several 
values of p. For p = 0, the errors on the fitted parameters are also given, and reflect the experimental uncertainties in 1171 . 
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-0.8 -0.6 



-0.4 -0.2 

p (GeV 2 ) 



0.2 



FIG. 1: x vs P 



K. 



(p=o) 



0.853 ±0.015 
0.278 ± 0.040 
0.394 ±0.020 
-0.125 ±0.021 
-0.157 ±0.023 



-0.011 ±0.024 
-0.497 ±0.009 
0.311 ±0.026 
0.550 ±0.024 
-0.595 ±0.017 



-0.476 ±0.018 
0.211 ±0.019 
0.817 ±0.025 

-0.247 ±0.055 
0.032 ±0.019 



-0.148 ±0.024 
0.774 ±0.013 

-0.102 ±0.022 
0.588 ±0.015 

-0.153 ±0.017 



0.156 ±0.017 
-0.181 ±0.035 
0.266 ±0.048 
0.525 ±0.020 
0.773 ±0.016 



(18) 



The fits are sensitive to the input mass of /o(1370) which has the largest experimental uncertainty in (|17|l . We 
find that very good % 2 fits can be obtained for 1.31GeV < m[/ (1370)] < 1.45GeV. Outside this interval, however, 
the goodness of the fits significantly reduces. This observation further restricts the m[/o(1370)] in i|17|) . For several 
values of m[/o(1370)] the fitted parameters are given in tabled and the corresponding rotation matrices are given in 
Appendix A. The resulting glueball mass is in the range 1.54 GeV < < 1.61 GeV. For this range of m[/o(1370)], the 
sensitivity of the glueball content of the /o(1370), /o(1500) and /o(1710) on the mass of /o(1370) are examined in figure 
121 and show a strong dependence. For example, we see that below the central value of m[/o(1370)] = 1.35 GeV, the 
glueball component of the /o(1710) dominates those of the /o(1500) and /o(1370). However, above to[/o(1370)] » 1.4 
GeV, the /o(1500) acquires the largest glue component followed by the components of the /o(1710) and /o(1370). 
The order again changes around m[/o(1370)] ~ 1.42 GeV. Of course this periodic-like behavior is not surprising as 
the rotation matrix, which depends on ten mixing angles, goes through periodic variations as we tune ro[/o(1370)]. 
Therefore, due to this sensitivity on the mass of /o(1370), a precises extraction of the glueball content of these states 
requires an accurate knowledge of the m[/ (1370)]. 

The correlation between the mixing parameter p and the input mass for the /o(1370) are also examined. Although 
the texture of the rotation matrix varies with these two parameters, the glueball mass remains more or less within 
the same intervals obtained by uncorrelated variations of p and m[/o(1370)] (see above). The overall numerical work 
shows 



1.47 GeV < m G < 1.64 GeV 



(19) 



in agreement with the lattice QCD estimates [23. There are values of p and m[/ (1370)], for which the glueball 
masses are as low as 1.32 GeV. However, for these cases the /o(980) acquires a large glueball component which is not 
consistent with either the molecule or the four-quark description of this state. 

It is important to note that the coupling of the glueball to both nonets N and N' should be taken into account. 
This is examined in table IIIII for p = and m[/o(1370)] = 1.35 GeV. The first column corresponds to the general case 



6 



0.8 



0.6 



0.4 



0.2 



A A 



1.32 1.34 1.36 1.38 1.4 
m[f (1370)] (GeV) 



1 .42 1 .44 



FIG. 2: The dependence on the m[/o(1370)] of the magnitude of the glueball component of /o(1370) (circles), /o(1500) (squares), 
and /o(1710) (triangles). 



Fitted 




m[/o(1370)] (GeV) 






Parameters 


1.34 


1.36 


1.38 


1.40 


c (GeV 2 ) 


0.143 


0.190 


0.193 


0.197 


d (GeV 2 ) 


-0.00692 


-0.00921 


-0.00968 


-0.0104 


c'(GeV 2 ) 


-0.00773 


-0.0161 


-0.0167 


-0.0215 


d! (GeV 2 ) 


-0.0104 


-0.00908 


-0.00699 


-0.00440 


e (GeV 2 ) 


-0.318 


-0.319 


-0.337 


-0.344 


/ (GeV 2 ) 


0.101 


0.125 


0.167 


0.188 


9 (GeV 2 ) 


1.291 


1.221 


1.204 


1.194 


ma (GeV) 


1.607 


1.563 


1.552 


1.546 



TABLE II: Best numerical values for the unknown parameters in Lagrangian 1131 . with p 
m[/ (1370)]. 



and several values of the 



where the coupling of the glueball to both nonets N and N' is taken into account. We see that scalar glueball coupling 
to N is larger than its coupling to N', In a recent work [2^ a similar analysis is given in which the coupling to the 
lowest-lying nonet is neglected. To illustrate the importance of the glueball coupling to both nonets, in the second 
and third columns of Table IIIII the glueball coupling to N and N' is respectively suppressed. We see that the result 
is sensitive to the glueball coupling to N and N'. For example, when e = the glueball mass significantly increases 
to 1.45 GeV, whereas when / = the glueball mass decreases to 1.22 GeV. (We see that the result seems to be more 
sensitive to e.) Other fitted parameters in Table ITTT1 also show a similar sensitivity to the suppression of the couplings 
e and /. Therefore, it is important to consider a general treatment in which both couplings are taken into account. 



IV. INTERACTION LAGRANGIAN 



To further examine the present model for the scalar mesons, we need to investigate the interaction Lagrangian for 
these states and study their partial decay widths to several two-pseudoscalar channels. We saw in previous section 
that the experimental input mass of /o(1370) is the main source of uncertainty on the rotation matrix which in 
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Fitted Parameters Fit 1 Fit 2 Fit 3 

c (GeV 2 ) 0.16 0.081 0.21 

d (GeV 2 ) -0.0078 -0.0067 -0.0098 

c'(GeV 2 ) -0.013 -0.027 -0.013 

d' (GeV 2 ) -0.0093 -0.0058 -0.013 

e (GeV 2 ) -0.32 -0.25 

/ (GeV 2 ) 0.12 0.13 

g (GeV 2 ) 1.27 1.45 1.22 



TABLE III: Numerical values for the unknown parameters in Lagrangian (1131) . The first column corresponds to the case where 
the scalar glueball couples to both the light scalar meson nonet V as well as to the heavy scalar meson nonet N . The second 
and third columns represent cases where the scalar glueball only couples to the nonet N' , and to the nonet V, respectively. 



turn affects the interaction Lagrangian and decay width calculations. Therefore, without an accurate knowledge of 
to[/o(1370)] we can only give a qualitative description of the decay properties. Moreover, although various decay 
ratios are recently reported by the WA102 collaboration [24[ (see table II . the experimental status of the individual 
decay widths is not quite clear. In PDG pj some of these decay widths are given, even though they are not used in 
any averaging. Nevertheless, in this section we present the interaction Lagrangian and give a preliminary study of 
several decay ratios. This will provide a basis in which the upcoming experimental data on the decay widths/ratios 
of the scalar mesons can be analyzed. 

The scalar-pseudoscalar-pseudoscalar interaction takes the general form: 



rl=0 



C 



7=1/2,1 



BTr (AO Tr 



-D'Ti (AO Tr (dpt) Tr (0„ 



) + DTt (AO Tr (0^0) Tr (0^0) + B'Tr (N') Tr (d^d^) 
EGTt (d^d^) + FGTr Tr 



(20) 



where B and D are unknown coupling constants describing the coupling of the four-quark nonet N to the pseudoscalars. 
Similarly, B' and D' are couplings of A' to the pseudoscalars. E and F describe the coupling of a scalar glueball to 
the pseudoscalar mesons, and C 1 ^^ 2 ' 1 is the interaction Lagrangian for / = 1/2,1 states given in © together with 
parameters in © . The pseudoscalar part of the Lagrangian can be found in ref . jj] • The interaction Lagrangian l|20l) 
can be rewritten as: 



where 7* s / is the coupling of the i-th. scalar [see Eq. (|15|l ] to pseudoscalars s and s' , and is given by 



(21) 



with K denned in (|15f) and 



E 

i 



hss>K) 1 



(22) 



„NS 



IS 



..INS 



Ts 



(23) 



in which the diagonal elements are the couplings of the pseudoscalars s and s' to the f^ s , /(f , /'q , f'o S [defined 
in (|14|l ] and the scalar glueball G, respectively. The diagonal elements for all decay channels ss' are listed in the 
Appendix B. 

To determine the unknown couplings we need to fit the prediction of this Lagrangian to experimental data. Here 
we use the estimates of the decay ratios by the WA102 collaboration |24[ in table HVl We should note, however, that 
the decay ratios alone are not sufficient to determine the free parameters and need to be supplemented by more data 
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such as the individual partial decay widths, or the total decay widths 

r tota/ [/ (980)] = 40—100 MeV 

r* ota '[/ (1370)] = 200 — 500 MeV m] 

r* ota '[/ (1500)] = 109±7 MeV 

r* ota '[/ (1710)] = 125 ±10 MeV 

For example, with p — and m[/o(1370)] = 1.35 GeV, a numerical study of the decay ratios is given in table 
IIVI together with the predicted decay widths in table E] The result is compared with estimates extracted from 
other works, and shows an overall qualitative agreement, even though some of the predicted decay widths such as 
r[cr(550) — > 7T7r] or r[/o(1500) — > tttt] do not quite agree with other investigations 1 . However we should again 
note that the current experimental status of the decay widths is not quite established, and that together with the 
large uncertainty on the mass of the /o(1370) are the main sources of error in estimates given in tables Hvl and Ivl 
Certainly more work is needed to investigate the parameters of the interaction Lagrangian in more detail, and study 
its correlation with the mass of /o(1370). We postpone this goal for future works. 



Fitted Parameters 



B (GeV -2 ) 


-1.0 - 


► -0.7 


B' (GeV" 2 ) 


0.8 - 


> 1.0 


D(GeV~ 2 ) 


-0.1 - 


-> 2.5 


D' (GeV 2 ) 


-3.5 - 


-> 0.5 


E (GeV" 2 ) 


-1.6 - 


> -0.3 


F (GeV" 2 ) 


-2.2 - 


-» 2.4 



Decay ratios This fit WA102 Collaboration 

1.4 — 6.5 2.17 ±0.9 

< 0.23 0.35 ±0.21 



r[/ (1370)^7T7r] 
T[f (1370)^KK] 

T[f (1370)^KK] 
r[f (1500)^KK] 

r[/ (1500)-> 

ri/otisoo)^^^] 
r[/ (i500)— 7777 

TSo^V °- 3 ^ L7 0.52 ±0.16 

r r g^ 0.19- 2.6 0.20 ±0.03 

r^S^t 0.06-0.21 0.48 ±0.14 

^'^r'"'' 1 0.04 — 0.16 < 0.05 



< 0.16 0.32 ±0.07 

> 5.3 5.5 ±0.84 



r\Mirio)-.r,r, 



TABLE IV: Numerical values for the parameters in the scalar-pseudoscalar-pseudoscalar Lagrangian [Eq. J3JJJ], obtained by 
fitting its prediction for decay ratios (with the specific choice of m[/o(1370)] = 1.35 GeV and p — 0) to the experimental data 
by WA102 collaboration. 



V. SUMMARY AND CONCLUSION 



In this work we studied the I = scalar mesons below 2 GeV [<r(550), / (980), / (1370), / (1500), and / (1710)] 
using a non-linear chiral Lagrangian which is constrained by the mass and the decay properties of the 1 = 1/2 and 
1=1 scalar meson below 2 GeV [k(900), Kq (1430), a (980) and an(1450)]. In this framework the lowest-lying four- 
quark scalar meson nonet N mixes with the next-to-lowest lying two-quark nonet N' and a scalar glueball G. We 
showed that this model can describe the mass spectrum of the scalars, and studied the correlation between the mass 
of / (1370) and the substructure of these states. We showed that consistency of this model with the experimental 
mass spectrum favors 1.31 GeV < m[/o(1370)] < 1.45 GeV, and sets a bound on the scalar glueball mass in the 
1.47 GeV to 1.64 GeV range. We also showed that it is important to take into account the coupling of the scalar 



Note that the F[a — > tttt], which is proportional to the 7 2 ?r7r , is not the same as the Breit-Wigner width that appears in the denominator 
of the <r propagator 1 1(1 . This deviation from a Breit-Wigner shape, is also a characteristic of the k(900) meson llll . 
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Decay Widths (MeV) This Model Extracted from other works 



T[a(550) - 


> 7T7r] 


53 -> 61 


w 100 [10] 


r[/ (980) - 


-> 7T7r] 


107 -> 116 


» 65 flOl 


r[/ (1370) 


— > 7T7r] 


121 -> 238 


34 -> 175 [25J 


r[/o(1370) 




36 -» 85 


44 -> 240 [25] 


r[/„(1370) 


-> 7777] 


< 20 




r[/ (i5oo) 


— > 7T7r] 


103 -> 329 


36 ^ 65 [25] 


r[/ o (i500) 


-> Kj?] 


2 -> 17 


2.4 7.5 [25] 


r[/ o (i500) 


-» ^ 


< 62 




r[/ (i500) 


-» W] 


< 23 


large [26J 


r[/ (i7io) 


— > 7T7r] 


< 92 


1.7 -> 5.5 [2JJ 


r[/ (i7io) 


-> A" A"] 


3 -> 36 


22 -> 64 [271 


r[/ (i7io) 


w] 


< 8 


6 -> 24 [27] 


r[/„(i7io) 


-> 7777'] 


< 1.5 





TABLE V: Partial decay widths of J = scalars predicted by the fit in Table llVl 



glucball to N as well as to TV'. We found that the <r(550) is mainly a non-strange four quark state, whereas the 
substructure of other 1 = states is sensitive to the mass of the Jo (1370). The numerical results show that the 
/o(1500) and /o(1710) have significant glueball admixtures. We also investigated the interaction Lagrangian and gave 
a preliminary study of the decay widths of the 1 = scalars into various pseudoscalar-pseudoscalar channels. Probing 
scalar-pseudoscalar-pseudoscalar couplings is important for low energy processes such as 77' — > 3n, 77' — » t]tttt, in which 
the scalar mesons are expected to play important roles |l2l l28j , and therefore are interesting directions for future 
works. It is also interesting to examine this model with higher order effects, such as higher derivative terms or more 
complex mixing terms between two and four quark nonets. 
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APPENDIX A: THE ROTATION MATRICES 

1. m[/ (1370)] = 1.35GeV and p a variable: 



K 



(p=-0.5GcV 2 ) 



0.829 


0.033 


0.030 


0.556 


0.034 


-0.289 


0.605 


0.648 


0.362 


-0.009 


-0.431 


-0.099 


-0.491 


0.692 


-0.292 


-0.066 


0.411 


-0.430 


0.049 


0.800 


0.197 


0.674 


-0.393 


-0.280 


-0.524 



(Al) 



"(p^-O^GeV 2 ) 



0.884 


-0.047 


-0.047 


0.463 


-0.003 


-0.224 


0.402 


0.685 


0.535 


-0.183 


-0.368 


-0.362 


-0.482 


0.616 


-0.349 


-0.128 


0.602 


-0.445 


0.266 


0.594 


0.129 


0.586 


-0.315 


-0.223 


-0.701 



(A2) 
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K 



(p=-0.2GcV 2 ) 



K 



(p=0.075GcV 2 



0.927 


-0.138 


-0.292 


0.155 


0.108 


0.039 


-0.107 


0.429 


0.783 


-0.436 


0.329 


0.342 


0.804 


-0.293 


0.210 


-0.141 


0.455 


-0.126 


0.524 


0.694 


0.104 


0.803 


-0.264 


-0.042 


-0.522 


0.842 


-0.016 


-0.490 


-0.183 


0.133 


0.274 


-0.531 


0.156 


0.769 


-0.167 


0.414 


0.062 


0.852 


-0.232 


0.212 


-0.077 


0.302 


-0.075 


0.435 


0.842 


0.199 


0.789 


0.057 


0.365 


-0.448 



(A3) 



(A4) 



2. p = and m[/ (1370)] a variable: 



K 



(m[/ (1370)]=1. 34GeV) 



K 



(m[/ (1370)]=1.36GeV) 



K 



(m[/ (1370)] = 1.38GeV) 



K 



(m[/ (1370)]=1.40GeV) 





0.844 


-0.003 


-0.490 


-0.159 


0.149 




0.298 


-0.506 


0.220 


0.762 


-0.163 




0.398 


0.348 


0.798 


-0.096 


0.274 




-0.129 


0.558 


-0.273 


0.602 


0.485 




_ -0.155 


-0.559 


0.028 


-0.148 


0.801 




0.853 


-0.034 


-0.476 


-0.142 


0.159 




0.277 


-0.467 


0.230 


0.784 


-0.196 




0.391 


0.267 


0.821 


-0.149 


0.283 




-0.126 


0.486 


-0.215 


0.554 


0.629 




0.167 


0.688 


-0.035 


0.192 


-0.679. 




0.867 


-0.030 


-0.447 


-0.124 


0.181 




0.241 


-0.481 


0.180 


0.785 


-0.251 




0.389 


0.204 


0.862 


-0.122 


0.220 




-0.110 


0.490 


-0.148 


0.570 


0.633 




0.167 


0.698 


-0.044 


0.170 


-0.674 _ 




0.883 


-0.031 


-0.417 


-0.100 


0.189 




0.197 


-0.483 


0.136 


0.793 


-0.283 




0.384 


0.132 


0.895 


-0.120 


0.138 




-0.085 


0.502 


-0.061 


0.567 


0.645 




0.163 


0.705 


-0.049 


0.158 


-0.671 



(A5) 



(A6) 



(A7) 



(A8) 



APPENDIX B: THE COUPLING CONSTANTS AND DECAY WIDTHS 

The coupling of I = states f^ s , Jq, f , J'q S and G to different two-pseudoscalar channels are: 



11 



T^r = (Bl) 

7 f/ = -V2(B-A) (B2) 

7 ix - -2(25 -A) (B3) 

-2^5 (B4) 



IKK 



7^ = -\/2(5 + I3) + -(C + 2A + 4£>)sin20 p - V2{C + D) cos 2 9 p (B5) 
1 



7^ S = -(B + D) + — (C + 2D)sm29 p -(A + D) cos 2 9 p - C sin 2 p (B6) 
V2 

-\/2(C + D) sin 29 p -{C + 2A + AD) cos 2(9 p (B7) 



lipi' 

y%j? = - v^(C + 2£>) cos 26» p - (A - C + D) sin 29 p (B8) 



-V2{B + D) - - (C + 2A + AD) sin 20 p - \/2{C + D) sin 2 9 p (B9) 

= _ (_g + _ ( C + 2j d) S i n 29 p ~ (A + D) sin 2 6> p - C cos 2 <9 p (BIO) 
V2 

7& - -V2(B'-A') (Bll) 

7^ S = -2B> (B12) 

T&r = -2V2S' (B13) 

= —2{2B' - A') (B14) 

7 'S = ~{B' + D') + —=(C' + 2D') sin 29 p - (A' + D') cos 2 P - C" sin 2 6» p (B15) 

= -v^(-B' + D') + ^(C' + 2A' + AD') sin 29 p - V2(C*' + D') cos 2 P (B16) 

= -V2(C + 2L>')cos26> p - (A' - C + D')sm29 p (B17) 

= -V2(C + D')sm29 p - (C + 2A' + AD')cos29 p (B18) 

V? = — (S' + £)') /=(C' + 2-D') sin20 p — (A' + D') sin 2 9 p — C' cos 2 9 p (B19) 

V2 



hp/ 



-V2{B' + D') - -{C + 2A' + AD') sin 29 p - y/2(C + D 1 ) sin 2 9 p (B20) 

-V2E (B21) 

-2V2£ (B22) 

-£-F(l+cos 2 6> p -\/2sin26g (B23) 

"fa = -F(sm29 p + 2V2cos29 p ) (B24) 

7$„» = +sin 2 6» p + \/2sin26' p ) (B25) 

where P is the pseudoscalar mixing angle defined as 



lr)'ri> 

I 7T7T 

IKK 
7 G 



7? \ = / cos p -sin 9 p \ ( {4>\ + 4>l)/V2 
rf I \ sin 9 P cos 9 P I \ cf>\ 



with 9 P fa 37° |2£ 



The two-body partial decay widths of physical states Fj are: 



m - tttt] = 3 I j [M 2 s - 2m 2 ] 2 (B27) 



2 

32ttM f , 



(B26) 



r[F, ^KK] = | J 7 ^f 2 ) [Af 2 t - 2m^J^ (B28) 
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r[F. 




vv] 



(B29) 



(B30) 



where q is the center of mass momentum of the final state mesons, and for a general two-body decay A — + BC is 
given as: 



[1] Review of Particle Physics, Phys. Rev. D 66 (2002). 

[2] D. Black, A.H. Fariborz, F. Sannino and J. Schechter, Phys. Rev. D 58, 054012 (1998). 

[3] E. van Beveren, T.A. Rijken, K. Metzger, C. Dullemond, G. Rupp and J.E. Ribeiro, Z. Phys. C 30, 615 (1986). 

[4] S. Ishida, M. Ishida, T. Ishida, K. Takamatsu and T. Tsuru, Prog. Theor. Phys. 98, 621 (1997). See also M. Ishida and 

S. Ishida, Talk given at 7th International Conference on Hadron Spectroscopy (Hadron 97), Upton, NY, 25-30 Aug. 1997, 

hep-ph/9712231 
[5] E.M. Aitala et al, Phys. Rev. Lett. 89, 121801 (2002). 
[6] R.L. Jaffe, Phys. Rev. D 15, 267 (1977). 
[7] J. Weinstein and N. Isgur, Phys. Rev. D 41, 2236 (1990). 
[8] N.A. Tornqvist, Z. Phys. C 68, 647 (1995). 

[9] V. Elias, A.H. Fariborz, Fang Shi and T.G. Steele, Nucl. Phys. A 633, 279 (1998); Fang Shi, T.G. Steele, V. Elias, K.B. 
Sprague, Ying Xue and A.H. Fariborz, Nucl. Phys. A 671, 416 (2000). 
[10] F. Sannino and J. Schechter, Phys. Rev. D 52, 96 (1995); M. Harada, F. Sannino and J. Schechter, Phys. Rev. D 54, 1991 

(1996); Phys . Rev. Lett. 78, 1603 (1997). 
[11] D. Black, A.H. Fariborz, F. Sannino and J. Schechter, Phys. Rev. D 59, 074026 (1999). 
[12] A.H. Fariborz and J. Schechter, Phys. Rev. D 60, 034002 (1999). 
[13] D. Black, A.H. Fariborz and J. Schechter, Phys. Rev. D 61, 074030 (2000). 
[14] D. Black, A.H. Fariborz and J. Schechter, Phys. Rev. D 61, 074001 (2000). 

[15] D. Black, A.H. Fariborz, S. Moussa, S. Nasri and J. Schechter, Phys. Rev. D 64, 014031 (2001). 
[16] N.N. Achasov and G.N. Shestakov, Phys. Rev. D 49, 5779 (1994). 
[17] J. A. Oiler, E. Oset and J.R. Pelaez, Phys. Rev. D 59, 074001 (1999). 
[18] S. Teige et al, Phys. Rev. D 59, 012001 (1999). 

[19] D. Black, M. Harada and J. Shechter, Phys. Rev. Lett. 88, 181603 (2002). 

[20] F.E. Close and N.A. Tornqvist, J. Phys. G 28, 249 (2002). 

[21] F.E. Close and A. Kirk, Phys. Lett. B 483, 345 (2000). 

[22] T. Teshima, I. Kitamura and N. Morisita, J. Phys. G 28, 1391 (2002). 

[23] N. Ishii, H. Suganuma and H. Matsufuru, Phys. Rev. D 66, 014507 (2002); 

Xi-Yan Fang, Ping Hui, Qi-Zhou Chen and D. Schutte, Phys. Rev. D 65, 114505 (2002); 

C.J. Morningstar and M. Peardon, Phys. Rev. D 60, 034509 (1999). 

J. Sexton, A. Vaccarino and D. Weingarten, Phy. Rev. Lett. 75, 4563 (1995); 

G. Bali et al., Phys. Lett. B 309, 378 (1993). 
[24] D. Barberis et al, Phys. Lett. B 479, 59 (2000). 

[25] D.V. Bugg, A.V. Sarantsev, B.S. Zou, Nucl. Phys. B 471, 59 (1996). 
[26] D.M. Aide et al., Phys. Lett. B 201, 160 (1988). 
[27] R.S. Longacre et al, Phys. Lett. B 177, 223 (1986). 

[28] A. Abdel-Rehim, D. Black, A.H. Fariborz and J. Schechter, Phys. Rev. D 67, 054001 (2003). 
[29] V. Mirelli and J. Schechter, Phys. Rev. D 15, 1361 (1977). 



\l\m\ - (m B + m c ) 2 ] \m\ - (m B - mpf] 
2mA 



(B31) 



